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ABSTRACT 

We present an investigation of potential signatures of the formation of multiple stellar populations 
in recently formed extragalactic star clusters. All of the Galactic globular clusters for which good 
samples of individual stellar abundances are available show evidence for multiple populations. This 
appears to require that multiple episodes of star formation and light element enrichment are the norm 
in the history of a globular cluster. We show that there are detectable observational signatures of 
multiple formation events in the unresolved spectra of massive, young extragalactic star clusters. We 
present the results of a pilot program to search for one of the cleanest signatures that we identify - 
the combined presence of emission lines from a very recently formed population and absorption lines 
from a somewhat older population. A possible example of such a system is identified in the Antennae 
galaxies. This source's spectrum shows evidence of two stellar populations with ages of 8 Myr and 
80 Myr. Further investigation shows that these populations are in fact physically separated, but only 
by a projected distance of 59 pc. We show that the clusters are consistent with being bound and 
discuss the possibility that their coalescence could result in a single globular cluster hosting multiple 
stellar populations. While not the prototypical system proposed by most theories of the formation of 
multiple populations in clusters, the detection of this system in a small sample is both encouraging and 
interesting. Our investigation suggests that expanded surveys of massive young star clusters should 
detect more clusters with such signatures. 
Subject headings: globular clusters: general - galaxies: star clusters - galaxies: stellar content 



1. INTRODUCTION 

Traditionally, globular clusters have been considered 
as simple stellar populations (SSPs) - produced from 
a single short burst of star formation, resulting in a 
population of stars with similar ages and abundances. 
This view follows naturally from the narrow color- 
magnitude diagrams observed in Galactic globular clus- 
ters. The well-defined main-sequence turnoffs limit any 
age spread, while the narrowness of features such as the 
red giant branch limit abunda nce spreads, at least in 
iron (e.g. lAshman fc Zepflll99cj and references therein). 
Further evidence that the spread in iron-peak abun- 
dances is small comes from detailed spectroscopic stud- 
ies of individual stars in Galactic globular clusters (e.g. 
A[Fe/H] <0.05 dex for stars in 19 Galactic clusters, 
iCarretta et al.l l2009aD . The few globular clusters that 
differ from this picture, by having clear spreads in iron- 
peak abundances among their stars (e.g. w Cen, M22, 
M54, Terzan5), are typically the most massive globular 
clusters and only represent a few percent of the over- 
all globular cluster popul ation (see e.g. the review's of 
iGratton et al. 2004, 20121 and references therein). 

However, it is now clear that multiple populations are 
the norm in the Galactic globular clusters. A key com- 
ponent in this revolution in understanding globular clus- 
ters is the significant star-to-star variation observed in 
light elements such as C, N, O, Na an d Mg (see e.g. 
ICarretta et alj|20lot IGratton et al.l [20121 and the many 
references therein) . Initial hints of variations in the abun- 
dances of light elements came from the spectra of small 
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samp les of giant stars observed in clusters (e.g. iKraftl 
1994). Subsequently, spectroscopy of much larger sam- 
ples of globular cluster stars, including main sequence 
stars, has found the same variations in light element 
abundances. The presence of these variations among un- 
evolved stars within the same globular cluster clearly in- 
dicates that these variations were present in the gas at 
the time the stars formed in the clusters. 

Notably, the variations in light elements such as C, N, 
O, Na and Mg have been seen within all of the many 
globular clusters now observed in suffici ent detail (e.g. 
ICarretta et alJl20Tol IGratton et al.ll2012l and the many 
references therein). Thus, these variations represent a 
fundamental aspect of globular cluster formation. More- 
over, the variations are not random but have clear pat- 
terns. Prominent among these is the well-known Na - O 
anticorrelation, such that stars with an overabundance 
of Na have an underabundance of O. A similar anti- 
correlation is seen in Al and Mg. Such anti-correlations 
are a signature of H-burning at high temperatures (e.g. 
ILanger et al.|[l9 93) , and thus place strong constraints on 
models for the origin of the material out of which the 
second population forms. It is also notable that stars 
with unusual patterns of light elements make up a sub- 
stantial fraction of all stars within a globular cluster. 
Specifically, about half of the stars in globular clusters 
have abundance ratios similar to those seen in halo stars, 
while the other half or so have different light element ra- 
tios, wi th for example a range of enhanced N and reduce d 
O (e.g. ICarretta et al.ll2009bl l2010t IGratton et al.ll20ll . 
Furthermore, there is a fairly smooth distribution of en- 
hanced Na and correspondingly depleted O among the 
stars with unusual abundances - from just slightly higher 
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[Na/O] than typical halo stars to an order of magnitude 
(or more) enhancement. 

A number of other signatures of multiple stellar popu- 
lations have been observed in the Galactic globular clus- 
ters. One example is the multiple or broadened main se- 
quences that are being obse rved in an increasi ng number 
of globular clusters (see e.g. lPiotto et al.ll2012l and refer- 
ences therein, for a discussion of clusters currently known 
to exhibit such signatures). These multiple sequences 
are often explained via a second generation of stars that 
have enhanced He abundance. Theories aimed at under- 
standing the ubiquitous light element variations observed 
in globular clusters generally predict such an enhance- 
ment. Spectroscopy of stars in the different sequences 
appears to confirm that they are associated with the ob - 
served abundance variations (e.e. iBragaglia et al.ll2010h . 
Variations in horizontal branch morphology have also 
been tied to multiple stell ar populations, par ticularly 
through variation in He (e.g. lGratton et al"1l2013l and ref- 
erences therein). Additional evidence for multiple popu- 
lations in star clusters has been presented based on HST 
obse rvations of Large M agellanic Clou d (LMC) clus- 
ters JM ackcv & Brobv Nielsen 2007: iMilone et al.l[200l 
iGoudfrooij et all 120111 iRubele et al.ll2013f ). These ob- 
servations demonstrated a broad or split main sequence 
turn-off in these intermediate age clusters, possibly sug- 
gesting an age spread of a few hundred Myr. 

The presence of multiple populations in clusters has 
challenged our understanding of their formation. The 
variation in light element abundances within every glob- 
ular cluster studied to date indicates that the mul- 
tiple population phenomena is a fundamental prop- 
erty of globular clusters. The nature of the varia- 
tions requires the second generation to be formed out 
of material enriched by hot H-burning. As discussed 
in the following section, the leading explanation ap- 
pears to be the formation of a second generation of 
stars out of ejecta from the cluster's asymptotic giant 
branch (AGB) sta rs mixed w ith so me pristine material 
fseee.g.lD'Antona fc yenturall2007t[D'Ercole et alJl2008t 
IConrov fc Gunnll2010HConrovll2012D ~ 

The goal of this paper is to help test ideas for the origin 
of the multiple populations by identifying and studying 
the signature of multiple populations in recently formed 
massive star clusters. Many such clusters have been dis- 
covered in the local universe, however most are too dis- 
tant to be resolved into individual stars. Therefore, our 
emphasis is to develop signatures of multiple stellar pop- 
ulations in young, unresolved star clusters. In Section^ 
we use stellar population synthesis models to investigate 
the signatures of multiple stellar populations on the inte- 
grated photometry (|2.1j) and spectroscopy (|2.2j) of a star 
cluster. Section l2~3l discusses the probability of detecting 
evidence of multiple stellar populations in observations 
of young, massive star clusters based on the main signa- 
ture we present. Finally, we present in Section [3J a pilot 
program to search for such signatures from star clusters 
in the Antennae galaxies. 

2. SIGNATURES OF MULTIPLE POPULATIONS IN 
UNRESOLVED STAR CLUSTERS 

The stellar populations of star clusters beyond our Lo- 
cal Group can not generally be resolved. Therefore, in 
order to study multiple stellar populations in extragalac- 



tic globular clusters, we need to identify their effect on 
the integrated emission from a cluster. In this section 
we consider the differences between the emission from a 
single stellar population and multiple stellar populations 
and discus their potential detectability. In the following 
discussion we take the simplest case of a multiple stellar 
population - that of a cluster with two distinct bursts 
of star formation. The age difference between the bursts 
of star formation in multiple population clusters is ob- 
servationally quite poorly constrained. Observations of 
star clusters in the Magellanic Clouds show spreads in 
the main sequence turnoffs that, if due solely to an age 
diffe rence, imply an age spread of several hund red Myr 
(e.g. IGoudfrooij et alJl20ll IRubele et~aII[20Tol) . How- 
ever, other effects such as rotation or He abundances may 
also act to broaden the main sequence. Therefore, this 
constraint represents an upper limit on the age difference 
between populations in these clusters. 

Theoretical considerations imply certain constraints on 
the age difference between the two populations. The 
spectroscopic evidence from the Galactic globular clus- 
ters requires that the second population must be formed 
from gas enriched by H-burning at high temperatures. 
Currently, the leading contender for producing the pol- 
luted second generation gas is intermediate mass (~4- 
8M Q ) AGB stars. In this scenario, hot H-burning occurs 
at the base of the convective envelope and is convected 
to the surface. This gas is then released into the in- 
tercluster medium via winds or with the ejection of the 
stars envelope. To produce the correct abundances in 
the second generation gas, this scenario requires AGB 
stars with masses in the range ~4-8M© and implies that 
the second ge neration should be ~40-200 Myr older than 
the first (e.g. iD'Antona fc Ventural 120071 ). Fast rotating 
massive stars (2O-12OAf , FRMS) have also been pro- 
posed as potential polluters. Here, hot H-burning occurs 
in the core of the star, its products are driven to the 
surface via rotational mixing and expelled via an equa- 
torial wind into the intercluster medium. However, this 
scenario, which invokes massive stars and hence short 
timescales, has difficulties in producing new stars which 
are not enriched in iron peak elements (from supernovae) 
and which have the correct light elem ent pattern (e.g. 
iRenzinil [20081 IConrov fc Spergelll201lD . Moreover, such 
models can not account for the possible age spreads sug- 
gested by the young LMC clusters. 

Spectroscopic observations of individual stars in the 
Galactic globular clusters have found that the second 
generation of stars comprise a large fraction of the to- 
tal stellar population in these clusters, around 50% (e.g. 
iCarretta et al.ll2009bf) . However, it should be noted that 
at earlier times the fraction of second generation stars 
is likely to have been much lower. Indeed, leading theo- 
ries for the formation o f multiple population s in clusters 
often requires this (e.g. lD'Ercole et al1l2008t ). This is be- 
cause the second generation gas is polluted by the first 
generation of stars, which makes it very hard to produce 
a second generati on which is >1 0% of the total cluster 
mass (e.g. iBekki fc Nor ris 2006). However, if the first 
generation of stars is preferentially lost, then this ratio 
can increase to that observed as the cluster evolves. Such 
a scenario may occur if the second generation of stars 
is more centrally concentrated than the first. Simula- 
tions have demonstrated that this process can produce 
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Fig. 1. — The evolution of a star cluster's U-B vs. V-I colors from the Bruzual & Chariot (2003) stellar population models. The black 
lines show a SSP with solar metallicity (aging from top left to bottom right). The blue dashed line shows a similar population, but with an 
additional burst of star formation, 100 Myr after the first. The second generation of stars is assumed to provide 10% of the cluster's total 
mass. The middle panel shows the same models but with different age spreads between the two bursts of star formation: 100 Myr (blue 
dashed line); 250 Myr (green dotted line); and 500 Myr (red dot-dashed line). The right panel shows SSP models with [Fe/H]=-0.33 dex 
(purple dot-dashed line), 0.00 dex (black solid line) and 0.56 dex (orange dotted line). Highlighted are the two regions in this color space 
where the multiple stellar populations have significantly different colors from SSPs. The arrow represents a reddening of E(B-V)=0.1 mag. 

with accurate photometry. Young, massive star clus- 
ters in the Antennae galaxies ar e found to lie in this 
region (e.g. lWhitmore et al.ll2010ft . However, several fac- 
tors could result in SSP clusters also having colors in 
this region. Firstly, reddening can cause SSP clusters of 
a few Myr to be shifted to these colors. Also, as can 
be seen from the right panel of Figure [I] SSP clusters 
with sub-solar metallicity can have similar colors. A 
broader multiwavelength approach may be able to ac- 
count for these f actors. However, there ar e also varia- 
tions between the iBruzual fe Chariot! (|2003t ) models and 
other commonly used st ellar population mode ls, such as 
th e Starburst99 m odel (|Leitherer et al.l I1999D and that 
of lMarastonl(|2005l ). While other SSP models predict the 
period of rapid reddening of the a cluster's V-I color, the 
tracks vary slightly and not all show SSPs avoiding re- 
gion 1. The other period where differences are observed 
between single and multiple population clusters is around 
a few hundred Myr (region 2 in Figure [lj. Here, massive 
giant stars from the second generation result in redder V- 
I colors than would be produced from a SSP cluster, but 
have little effect on the cluster's U-B color. A different 
metallicity SSP cluster could not explain such a shift, 
although reddening of a SSP cluster potentially could. 
Such a shift could be reliably measured, but only if the 
age difference between the two stellar populations was 
large (even at 500 Myr, the shift is only VT=0.2 mag). 
Interestingly, no young massive star clusters (with My<- 
10) in the Antennae or M83 ar c found with such colors 
(jWhitmore et al.ll20loT: IChandar et all 12010ft . A dearth 
of clusters in this region may argue against multiple stel- 
lar population clusters with large age spreads of more 
than a few hundred Myr (in agreement with current for- 
mation theories). 

We conclude that the identification of clusters with 
colors in regions 1 or 2 may be suggestive that they 
host multiple stellar populations. However, complica- 
tions due to reddening and model/photometric accuracy 
mean that this combination of U-B and V-I colors on 



a cluster with similar fractions of fi rst to second genera- 
tion stars a fter a Hubble time (e.g. iD'Ercole et alJl2008t 
lBekkill20ll . 

2.1. Broadband optical colors: U-B vs. V-I 

In Figured] we plot the broadband U-B and V-I colors 
of a st ar cluster, produced using the IBruzual fc Charlotl 
(2003) stellar population synthesis code. In the left 
panel, we consider the tracks traced by a solar metal- 
licity stellar population as a function of its age. The 
black line shows a SSP formed from a single burst of star 
formation. The cluster colors trace a path from top left 
to bottom right as the stellar population evolves. The 
biggest departure from this general trend occurs after a 
few Myr, when the cluster reddens rapidly in V-I due 
to the evolution of massive stars to supergiants. These 
models are found to be in good agreement with the ob- 
served U-B and V-I col ors of young massive st ar clusters 
in both the Antenna e (jWhitmore et all 12010ft and M83 
(|Chandar etal.ll2OT0ft . 

The blue dashed line in Figure [T] shows the color of 
the same cluster, but with an additional burst of star 
formation 100 Myr after the first. For these illustrative 
purposes we consider this second population to be formed 
with 10% of the total cluster mass, this is in the range 
predicted at early stages in the cluster's evolution (as dis- 
cussed above). The second panel shows the same cluster, 
but with different times between the bursts of star for- 
mation (100, 250 and 500 Myr). It can be seen that by 
ages >1 Gyr, the two populations in these clusters sum 
to give very similar U-B and V-I colors to a cluster with 
a SSP. Therefore, such observations at these ages would 
be unable to distinguish a multiple from a single stellar 
population. 

Some differences can be observed in this color space. 
For a few Myr after the second burst of star formation 
the cluster colors become bluer and fill region 1 of the 
color magnitude diagram. Such clusters would lie off 
the tracks of a SSP and this shift would be detectable 
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Fi g. 2. — Spectra of a s olar m etallicity star cluster produced from 
the Bruzual & Chariot (2003) models. The red points show the 
spectrum of a cluster with only one burst of star formation. The 
black lines show the spectra of a similar cluster but but with an 
additional burst of star formation 100 Myr after the initial burst. 
Similarly to the cluster colors plotted in Figure 1, the second gen- 
eration is taken to be 10% of the total cluster mass. The cluster 
is shown at 100 Myr (when only one stellar population exists), 
250 Myr, 500 Myr and 1 Gyr (from top to bottom, respectively). 

their own can not conclusively demonstrate the presence 
of multiple populations in a cluster. A better under- 
standing of the variations between different stellar popu- 
lation models and the addition of photometry extending 
deeper in to the ultraviolet and infrared wavelengths may 
help to draw more conclusive evidence based on broad- 
band photometry. While an interesting and important 
prospect for the future, such work is beyond the scope of 
the current paper. Finally, we note a caveat to these color 
considerations is that the second population is proposed 
to be He enhanced. Such an enhancement may influ- 
ence the integrated cluster colors, by making the second 
generation stars bluer. Unfortunately, stellar population 
models which include He enhancement are still in their 
infancy and not y et complete enough to cons ider in a sim- 
ilar fashion. The iBruzual fc Chariot! (|2003|) models also 
exclude the contribution of nebular emission to the clus- 
ter's broadband colors. Investigation of the Starburst99 
mode ls (which include nebular emission, iLeitherer et al.l 
19995 suggests that such emission can influence the clus- 
ter colors slightly at very young ages but should not sig- 
nificantly effect the general signatures of multiple popu- 
lations discussed. 
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Fig. 3. — Spectra of a SSP star cluster with solar metallicity 
shown at 4, 8, 12 and 20 Myr (from top to bottom). These spectra 
are produced using the Starburst99 high resolution models with the 
addition of nebula emission lines (see text for details). The spec- 
trum of the 4 Myr cluster is dominated by strong nebular emission 
lines. At 8 Myr these lines are weaker, but still clearly detectable. 
However, by 20 Myr the evolving stellar population produces very 
few ionizing photons and the spectrum is dominated by the stellar 
continuum, featuring strong absorption lines. 

2.2. Optical spectroscopy 

Integrated spectroscopy of star clusters is available/ 
obtainable for many nearby galaxies. In Figure [2J we 
show the predicted spectra of a st ar cluster produced 
from t he stellar population models of lBruzual fc Charlotl 
( 2003) . The red circle points show the spectra of a solar 
metallicity SSP at four different ages - 100 Myr, 250 Myr, 
500 Myr and 1 Gyr. Overlaid are the spectra of the same 
star cluster at these times, but with an additional burst 
of star formation at 100 Myr. Again, we assume that 
this second population of stars comprises 10% of the to- 
tal mass of the cluster. In the 250 Myr old cluster, the 
multiple population cluster is slightly brighter, particu- 
larly at short wavelengths. However, given typical errors 
on observed spectra at these wavelengths and model de- 
pendent variations, such a spectrum is adequately fit to 
a slightly lower mass and younger single stellar popu- 
lation model. Therefore, as observed in the colors pre- 
sented in Section 12. 1[ we find that the spectra of a mul- 
tiple stellar population cluster at these wavelengths and 
ages can be reasonably well represented as a SSP. More 
detailed consideration of the relative absorption lines 
strengths, rather than just considering the entire spec- 
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Fig. 4. — Starburst99 models with emission lines for solar metallicity SSPs at 4 (blue, left), 10 (blue, right) and 100 Myr (red). The lower 
panels show the spectra expected from a cluster composed of a 100 Myr stellar population and an additional minority population of 4 (left) 
and 10 Myr (right). We consider the initial mass ratio of the second population to initial population be 3, 5 and 10% (middle to bottom, 
respectively). In all cases, the absorption features of the old population and the emission lines from young population can be observed. 



trum with no weighting, may allow multiple populations 
to be distinguished (s uch as the methods used in inves- 
tigating galax ies, e.g. IVergelv et af] 2002; P anter et al.l 
I2003t IMathis et all [20061 ). These methods require high 
signal to noise spectra, reliable stellar population mod- 
els and a detailed understanding model limitations and 
fitting procedures. Such considerations are beyond the 
scope of this paper. Instead, we focus on a clearer signa- 
ture of multiple populations, discussed below. 

At ages of a few Myr the emission from a cluster fea- 
tures both nebular and stellar emission. This is because 
the young massive stars in the cluster can photoionize 
surrounding gas, resulting in a nebula that significantly 
influences a cluster's emission. Man y of the commonly 
used stellar populatio n models (e.g. iBruzual fc C hariot 
2003; Maraston 2005) only consider the stellar emis- 
sion from a cluster. Currently, only two models include 
bot h stellar and nebula r emission. The GALEV mod- 
els (|Kotulla et al.|[20 09h are the only models which cur- 
rently predict spectra for a stellar population with neb- 
ular continuum and emission lines included. However, 
these models focused on the effects of nebular emission on 



photometric colors and currently produce spectra with a 
resolution of only ^20A, lower than most observed spec- 
tra. 

The other model to consider the ionization of gas by 
a stellar population is the Starburst99 model. This code 
predicts the colors and spectra for a given mass, age and 
metallicity. While high resolution spectra (~2.5A) can 
be produced by this code, they only include the contin- 
uum contribution from any ionized gas, excluding any 
emission lines that would be produced. However, Star- 
burst99 does compute the luminosity of the H^ emis- 
sion line. We use this to predict the luminosity of all of 
the Balmer lines using the relative line strengths quoted 
by lOsterbrockl (|1989t assuming Case B, a gas density 
of 10 2 cm -3 and T=10 4 K). Additionally, we estimated 
the strength of the other prominent lines expected, the 
[O n] and [O in] lines. The combined strength of these 
oxygen lines was estimated by assuming solar metallicity 
and using the empirical relationship between metallicity 
and R 23 = (TO ill A 3727 +[0 nilAA 4959,5007 )/R , as 
presented in lKewlev fc Dopital ((2002). The ratio of [O n] 
to [O in] was estimated as [O n] =2.5[0 in] , assuming 
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a typical ionization parameter (Ke wlev fc Dop ita 2002). 
Finally, we fixed [O in] A 5007 = 3x[0 in] A 4959. This 
last relation is required by atomic physics, but it should 
be noted that some of the other ratios can vary. For 
our illustrative purposed, we have simply used typical 
values for these ratios. The ionized gas in HII regions 
generally has a low velocity dispersion and would pro- 
duce unresolved lines in spectra with the resolution of 
these models. We therefore added these emission lines to 
the Starburst99 synthetic spectra as Gaussian functions 
with FWHM=2.5A, similar to the spectral resolution of 
the models. 

Figure [3] shows the high resolution Starburst99 model 
spectra for a solar metallicity SSP with emission lines 
included. We show the model at four ages (4, 8, 12 and 
20 Myr). Massive, young stars in the 4 Myr stellar pop- 
ulation produce many ionizing photons and the cluster's 
spectrum features strong emission lines. By 12 Myr the 
nebular lines are much weaker due to the evolution of the 
most massive, hottest stars and subsequent reduction in 
the number of ionizing photons. However, ionizing pho- 
tons are still present and the stronger lines are still clearly 
observable above the cluster's continuum. By 20 Myr the 
nebular emission is extremely weak and the spectrum is 
dominated by stellar emission. 

In Figure U we use these models to show the spectrum 
produced by adding a minority second stellar popula- 
tion of 4 Myr (left, blue) and 10 Myr (right, blue) to a 
100 Myr old population (red). The initial ratio of first 
to second generation stars at these young ages is still 
poorly understood. We therefore consider three different 
cases where the minority (younger) population comprises 
3, 5 and 10% of the cluster's mass. In all cases the stellar 
like spectrum from the older stellar population, featuring 
broad absorption lines, can be observed in the spectra of 
the multiple population clusters. In addition to this, the 
narrow emission lines produced by the younger second 
population can also be observed. For the 4 Myr case, 
the emission lines are extremely strong and all of them 
can be observed in the multiple population cluster, even 
when the younger population comprises only 3% of the 
total population. When the second population reaches 
10 Myr, the nebular emission is much weaker. However, 
the strongest emission line (H Q ) is still observable above 
the stellar like continuum. Therefore, assuming that the 
second population in a two population cluster is younger 
than ~10 Myr and comprises more than a few % of the 
total cluster population, the spectrum of the multiple 
population cluster should feature both strong absorption 
and emission lines. 

Could a SSP produce such a spectrum? The ratio of 
emission to absorption line strengths in the spectrum of 
a young cluster (<10 Myr) would be too high to fit such 
a spectrum. Similarly, the emission lines in the spectrum 
of the multiple population cluster could not be produced 
by a stellar population older than ^20 Myr. We note 
that older populations can potentially feature emission 
lines from planetary nebulae. However, t hese are known 
to be extremely rare in globular clusters (Peacock ct al. 
12012ft and would struggle to produce both the predicted 
emission line flux and the correct ratio of H Q /[O in] . 
There is a small range of ages around 12-16 Myr where 
a SSP cluster may feature absorption and emission lines. 
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Fig. 5. — Starburst99 spectra of a SSP cluster at 14 Myr (blue 
dotted line) and a multiple stellar population cluster with 90% of 
its stars at 100 Myr and 10% at 10 Myr (red solid line). The full 
cluster spectrum is shown in the top panel and lower panels show 
zoomed in sections. The two models have similar absorption and 
emission lines. However the continuum of the multiple population 
is dominated by the more numerous older population of stars and 
is much redder than the 14 Myr cluster. 

Such a cluster could not fit the strong emission lines in 
the spectrum of a multiple stellar population cluster with 
a second generation younger than ^8 Myr. However, it 
could have a similar ratio of absorption to emission lines 
as a multiple stellar population with an older second gen- 
eration (~10 Myr). In Figure [5] we compare two such 
models (scaled to the same flux at 5000A) - a single 
stellar population at 14 Myr and a multiple stellar pop- 
ulation with 90% of the cluster mass in stars that are 
100 Myr old and the other 10% in stars that are 10 Myr 
old. It can be seen that the two models can produce 
similar strength emission and absorption lines. However, 
the continuum of the multiple stellar population (dom- 
inated by the 100 Myr population) could not be well 
represented by the younger 14 Myr population - with 
the 14 Myr population producing too much flux at short 
wavelengths and too little flux at long wavelengths. This 
is analogous to the gap in region 1 of the color-color plot 
shown in Figure [TJ Such a discrepancy can be used to 
distinguish between these intermediate age, single stel- 
lar population clusters, and a multiple stellar population 
cluster. We conclude that, unlike at older ages, the spec- 
trum of a multiple stellar population cluster, where the 
younger population is <10 Myr, can not be well repre- 
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sented by a single stellar population cluster with either 
a young, old or intermediate age. Consideration of these 
models therefore predicts a key period when a cluster 
comprised of multiple stellar populations can be clearly 
identified based on its integrated spectroscopy: 

When the second stellar population in a star cluster 
is younger than ~10 Myr, the presence of multiple stel- 
lar populations should be identifiable in integrated spec- 
troscopy due to a stellar continuum that indicates an 
older stellar population in addition to emission lines, that 
require an additional, very young population. 

When the younger population is older than ~10 Myr, 
any nebular emission is very faint and the key spectral 
features of the two populations are similar. This makes 
clear identification of more than one stellar population in 
the observed integrated emission much more challenging 
and more model dependent. 

2.3. Probability of observing such signatures 

Given the above considerations, we believe that the 
clearest detectable signature of multiple stellar popula- 
tions in an unresolved cluster occurs when the youngest 
population is <10 Myr. The presence of emission lines 
from this population in addition to the older cluster spec- 
trum should be detectable in spectroscopy of the cluster. 
This requires an age spread of >20 Myr between the 
populations (so that the first population is dominated 
by absorption features). This would be expected if the 
second generation gas is enriched via winds from AGB 
stars, as has been proposed. 

We can estimate the likelihood of detecting a cluster 
during the period when the multiple stellar populations 
can be detected using this method. To do this, we assume 
that all clusters have two stellar populations (additional 
bursts of star formation will tend to increase the like- 
lihood of detection by increasing the number of epochs 
when the young population produces nebular emission). 
The time between the two bursts of star formation (and 
whether this is universal to all clusters) is currently un- 
known. We therefore assume that the second population 
has an equal probability of being formed at any epoch 
in the range permissible by the leading multiple popula- 
tion formationary theory, that invoking enrichment from 
AGB stars - 40-200 Myr. Finally, using this method, 
the second generation is detected via its associated neb- 
ular emission and is therefore detectable for a period of 
~10 Myr. Given these considerations, we predict that 
there is a 6.5% chance of observing a random star clus- 
ter (with an age of 40-200 Myr) during the period when 
its second stellar population is younger than 10 Myr and 
thus detectable. Therefore, one would expect 1 in every 
15 such clusters to exhibit such a signature and there 
would a >99% chance of detecting such a signature from 
at least one cluster in a sample of 70 clusters. Spectro- 
scopic samples of this size are potentially observable in 
nearby galaxies. Candidates for clusters that are cur- 
rently at this interesting period in their evolution could 
be identified in photometric surveys as having colors sim- 
ilar to region 1 in Figure [T] 

3. OBSERVING MULTIPLE POPULATIONS IN MASSIVE 
YOUNG CLUSTERS: A PILOT STUDY 

3.1. Young massive clusters in the Antennae galaxies 



The interacting Antennae galaxies provide a rich 
source of young star clusters. These clusters have been 
the focus of many photometric and spectroscopic studies. 
In this section, we discuss a search for such signatures of 
multiple stellar populations in published spectroscopy of 
these clusters. 

A spectroscopic survey of 15 of the massive young 
star clusters in the A ntennae galaxies was presented by 
iBastian et al.l (|2009[ ). This survey confirmed that these 
clusters have ages in the range of a few Myr to a few hun- 
dred Myr and hence lie in the range where signatures of 
of multiple populations can be detected. We briefly re- 
view the data used by this study, but for a more d etaile d 
discussion we refer the reader to IBastian et al.l (2009) . 
The spectra of these clusters were obtained using the 
multi-object mode of the Gemini Multi Object Spectro- 
graph (GMOS) on Gemini-North under the program GN- 
2003A-Q-33 (PI: Gelys Trancho). Clusters were targeted 
based on preimaging supplemented by HST observations 
to confirm their cluster status. Masks with slit widths of 
0.75" were used to observe the clusters using the B600 
grating, with a resolution R=1688 at 461 0A. These spec - 
tra were presented and discussed in Bas tian et al.l (2009) 
and we do not repeat their analysis. However, we note 
that seven of the clusters targeted show emission line 
spectra consistent with young SSPs of a few Myr and 
seven clusters have spectra that are dominated by ab- 
sorption lines a nd consistent with SS Ps of a few 100 Myr 
(see figure 3 of IBastian et al.ll2009f) . Interestingly, the 
other cluster targeted (Till) shows both absorption and 
emission lines. While the absorption features can be well 
fit by a SSP of 80 Myr, the strong emission lines can not 
be produced by such a population. This source there- 
fore exhibits the signature expected from multiple stellar 
populations. 

3.2. The Antennae star cluster Till 

To further investigate the intriguing cluster Till, 
we obtained the raw spectra from the Gemini Science 
ArchivqJ. The spectra were reduced using the standard 
iraf Gemini packages. The raw spectra were first pre- 
pared for the Gemini packages using GPREPARE and nor- 
malized flatfield images were produced using GSFLAT. 
The edges of each slit on the image were found from the 
combined flatfield images using GSCUT. The task GSRE- 
duce was then used to subtract off the bias image (pro- 
vided by the Gemini archive), clean for cosmic ray hits, 
mosaic the three detectors together, apply the flatfield 
correction and cut the image (containing all of the multi- 
object spectra) in to individual spectra. A wavelength 
calibration was applied to the images based on associ- 
ated CuAr arc images, using the task GSWAVELENGTH. 
This produced the final set of six 2D spectral images of 
Till. The spectra were then located, traced and ex- 
tracted from these images using the iraf/noao APALL 
task with an extraction width of 9 pixels (~1.3"). Be- 
fore combining, the resulting spectra were then corrected 
for slit losses due to atmospheric differential refraction. 
This was required because the spectra were generally not 
taken at the parallactic angle - although the choice of the 
(fixed) slit angle minimized t his effect. This co rrection 
was based on the method of iFilippenkd (|1982f ) for our 

1 http://www3.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/gsa/ 
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Fig. 6.— CMOS spectrum of the cluster Till (black line, top). 
The middle panel compares this to an 80 Myr, solar metallicity SSP 
model from the Starburst99 code (red line). The bottom panel 
shows the residual after subtraction of the 80 Myr model. The 
model provides a good representation of the absorption line spec- 
trum. However, strong emission lines that are inconsistent with 
such a model are also present. Marked in this lower panel are the 
(rcdshiftcd) wavelengths of the Balmer series, ionized oxygen [O II] 
and [O III] , and neutral He [He i] lines. 

0.75" x 7.0" slit, measured airmasses of 1.3-1.5, seeing 
of 0.8-1.0"and angles to the parellactic angle of 3.5, 5.0, 
6.6, 18.5, 24.0, 43.1°. The maximum correction to the 
flux at 3700A, relative to 5000A, was a factor 1.5. 

The spectrum of Till is known to be well represented 
by a SSP model of so lar metallicity and an age of 80 Myr 
(jBastian et al.l 2009). We use such a model to apply an 
approximate flux calibration to the spectrum. This was 
done by fitting a smooth function to the ratio of the 
model to the GMOS spectrum and scaling the spectrum 
using this function. The final spectrum of Till is shown 
in the top panel of Figure [6) The s pectrum is in e x cellen t 
agreement with that presented by iBastian et al.l ([2009). 
The middle panel of this figure compares this spectrum 
with a Starburst99 SSP model with solar metallicity 
and an age of 80 Myr. The bottom panel of Figure 
[5] (blue) shows the residual emission after subtracting 
this model from the spectrum. While the model pro- 
vides a good representation of the cluster's spectrum, 
strong emission lines are also present. These are nebu- 
lar emission lines with the Balmer series, ionized oxygen 
O III A 3727 and [O m]AA 4959,5007 , and neutral He 
He i] AA4471,5876 lines all clearly visible. The observed 
wavelengths of these lines are in excellent agreement with 
the expected rest frame wavelengths, assuming a velocity 
of 1383 kms -1 (as marked in Figure [5]). The luminosity 
of these emission lines was measured using SPLOT. They 
are well represented by a Gaussian with FWHM=3.5A, 
consistent with the spectral resolution of the data, as 




Fig. 7.— HST WFC3/UVIS images of Till in the broadband 
F336W (U), F555W (V) and F814W (I) filters. The bottom right 
panel shows the ACS/WFC narrowband F658N (H a ) image of this 
region and the location of the GMOS slit. The bottom left figure 
indicates the orientation of all of these images and highlights the 
two clusters, Tlll-1 and Tlll-2. 

expected for a typical HII region. 

3.2.1. The two components of Till: a binary cluster? 

Multiwavelength archival HST photometry of this clus- 
ter is available from several surveys of the Antennae. We 
utilized the most recent optical survey of the Antennae 
(obtained under PID 11577, P.I. Whitmore). This pro- 
vides broadband F336W, F555W, F625W and F814W 
photometry (from WFC3/UVIS observations), broad- 
band F435W and narrowband F658N photometry (from 
ACS/WFC observations). We obtained the pipeline re- 
duced and combined images of the cluster in these dif- 
ferent filters from the Hubble Legacy Archive. In these 
images, the cluster was found to be offs et by 2" from the 
location quoted by IBastian et al.l (|2009[ ). lying in our im- 
ages at a=12:01:53.3, ^=-18:51:37.8. This is likely due to 
an error in the absolute astrometry of either our or their 
images. However, the location we take for the cluster 
was confirmed to be the object targeted by the Gemini 
spectrosco py based on the image presented in figure 4 of 
IBastian etaJ] (|2009h . 

The narrowband F658N HST photometry (covering 
H Q at the velocity of this cluster) showed that Till 
does not exhibit strong excess emission in the H a fil- 
ter, as would be expected given the strong Balmer emis- 
sion lines observed in its spectrum. However, a source 
very close to the targeted cluster was found to have 
a strong H Q excess. At only 0.62" from the cluster, 
this additional source would be blended with the clus- 
ter emission in the Gemini observations and included in 
the 0.75" slit (see Figure [7]). We believe it to be ex- 
tremely likely that this second source provides the opti- 
cal emission lines observed in Tlll's spectrum. Both of 
the sources are resolved in these images and are believed 
to be two distinct clusters with absolute V-band mag- 
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nitudes of - 10.6±0.1 and -8 . 6±0.2 (assuming a distance 
of 22 Mpc, iSchweizer et al.lfeOOSl ). The photometry of 
these two clusters was performed using iraf/apphot. 
For the brighter, older cluster (Tlll-1) we used an aper- 
ture with a radius of 0.3". For the fainter, younger clus- 
ter (Tlll-2), we used a smaller 0.2"aperture because of 
contamination from nearby sources. Comparing the col- 
or s obtained for these two c lusters with the SSP models 
of iBruzual fc Charlotl ((20031) . we find that Tlll-1 is well 
represented by an 80 Myr old, solar metallicity cluster 
with mass 4O±4xlO 4 M and Tlll-2 by an 8 Myr old, 
solar metallicity cluster with mass 10±2xl0 3 M©. 

The HST photometry of the Till system therefore 
suggests that the old and young components observed 
in the GMOS spectrum are likely due to two distinct 
clusters. We conclude that Till is not an example 
of a single young cluster hosting multiple populations. 
However, the proximity of these two clusters to each 
other (0.62", corresponding to a projected separation of 
59 pc at the distance of the Antennae galaxies) raises 
the question of whether their coalescence could result 
in a single cluster with multiple populations. Binary 
clusters ha ve been proposed in other galaxies, such as 
the LMC (IDieball et alj [200l IMucciarelli et all [20121) 
and NGC 5128 (|Minniti et all 120041 ) and models show 
that binary clusters can be formed during high velocity 
cloud interactions, a s would be expected in th e Antennae 
(jFuiimoto fc Kumail Il997t IBekki et afl 12001 . The coa- 
lescence of star clusters potentially provides a method of 
producing a single globular cluster hosting multiple stel- 
lar populations - although such a scenario struggles to 
produce the observed abundance patterns (such as the 
similar [Fe/H] and the Na-0 anticorrelation) . 

The close alignment of the Balmer emission and ab- 
sorption lines in Figure [6] demonstrates that Tlll-1 and 
Tlll-2 must have similar radial velocities. Using the 
brightest emission lines in the spectrum (H,g , H 7 , H,5, 
[O m]AA 4959,5007 and [O n]A 3727 ) we measure Tlll- 
2's radial velocity, v = 1383 ± 4 kms -1 (where the error 
is estimated as the standard deviation of the velocities 
measured for each line). The radial velocity of Tlll-1 
was measured from the absorption lines using the cross- 
correlation iraf routine FXCOR. For our reference spec- 
trum we used the high resolution Starburst99 SSP model 
for a solar metallicity, 80 Myr old star cluster. The 
cross correlation was run over the regions of the spec- 
trum that were unaffected by the emission lines. We 
experimented with several different wavelength ranges, 
achieving consistent results. The most accurate correla- 
tion was produced by using most of the spectrum, from 
3900-6100A(excluding the emission lines), for which the 
cross correlation resulted in a radial velocity and asso- 
ciated error of 1391 ± 17 kms -1 . Because the stellar 
population models vary slightly in their resulting spec- 
tra, we also compared our spectrum to that o f a similar 
cluste r produced using the u pdated models of lMarastonl 
(|2005D ; [Maraston fc Strombackl (|2011[ ). Excellent agree- 
ment was found between the two models with the cross 
correlations producing velocity differences of only a few 
kms -1 . 

For Tlll-1 and Tlll-2 to be gravitationally bound, 
we require that the potential energy is greater than the 
kinetic energy, or: 



< 



R 



(1) 



Here M sys is the total mass of the two clusters, R is 
their separation and v is their relative velocity. Unfor- 
tunately, it is only possible to measure the line of sight 
velocity (vi os ) and the projected separation (R p ) between 
the two clusters. We substitute thes e parameter s in to 
Equation[T]by following the method of lDavis et al.l (1995) 
and defining the system to be at an angle a to the sky: 



V los R P 



2GM. 



< sin 2 acosa (0 < a < 90°) 



(2) 



*</■-• 



For the Till system, we have measured Mi=10 M , 
M2=1O 4O M and i? p =59 pc, but the angle a is unknown. 
However, for most angles (10 < a < 80°), we find that 
the maximum permissible vi os between two bound clus- 
ters is 3.1 < vi os . esc < 4.8 kms -1 (depending on a). The 
measured line of sight velocity difference between Tlll-1 
and Tlll-2 is consistent with the clusters being bound. 
However, the accuracy of the radial velocity estimates is 
too low to be conclusive. 

3.2.2. Abundance Analysis of Tlll-1 and Tlll-2 

The spectrum of Tlll-1 suggests that the metallic- 
ity of its stellar pop ulation is around solar, 0.9±0.3Z Q 
(jBastian et al.l 120091 ) . The metallicity of the younger 
cluster, Tlll-2, can also be estimated based on its emis- 
sion line strengths. To do this we use the R23 diagnostic, 
defined as: 

[O n]AA 3726, 3729 + [O m]AA 4959, 5007 , , 

«-23 = 7Z (3) 

From the continuum subtracted spectrum (presented 
in the bottom panel of Figure EJ, we estimate these 
line fluxes as: Up =4.6xl0 37 ergs" 1 ; [O m]AA 4959,5007 
^^xlO^ergs- 1 ; [O 11] AA 3726,3729 =5.0xl0 37 erg S - 1 . 
This implies log(R 23 )=0.25 for Tlll-2. Based on the 
noise of the continuum, we estimate that the formal er- 
rors on these strong emission lines are <1%. The uncer- 
tainties are likely to be larger than this due to the flux 
calibration of the spectrum, particularly for the ratio of 
[O 11] to the redder lines. However, the uncertainty on 
the line fluxes are believed to be smaller than the large 
systematic error associated with using these lines to find 
[Fe/H] and the He abundance, discussed below. 

Empirical relationships between R23 and log(0/H) 
have been extensively studied, although significant scat- 
ter is known t o exist producing sy stematic errors 
of 0.2-0.5 dex (jKobulnickv fc Kewlevl [200l and ref- 
erences therein). We use the relationship between 
R23 and log(Q/H) that is p resented in equation 16 of 
iKobulnickv fc Kewlevl (|2004D . This is the average of 
the hi gh me tallicity relationships pres ented bv lMcGaughl 
(119911 ) and iKewlev fc Dopital ((2002). Using this rela- 
tion we find that log(O/H)+12=9.0 which, assuming so- 
lar abundances, implies a metallicity of 1.2Zq. This is 
higher than that of Tlll-1, but consistent given the pre- 
dicted uncertainties. 

In addition to the [O 11] , [O ill] , and Balmer emis- 
sion lines, [He i]A 4471 and [He i]A 5876 emission lines 
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are clearly seen in the spectrum (Figure 15]). In principle, 
these lines can be us ed to constrain t h e He a bundance in 
the Til 1-2 nebula. iBeniamin et al.l (|1999T) present the 
predicted emission coefficient, Airfu ne /n e nii e , as a func- 
tion of the gas density and temperature. Based on these 
predictions (for n e =10 2 cm -3 and T=10 4 K) and using 
the observed emission line flux ratio, /He, 4472//H fi =0.075, 
we can estimate the He number abundance as: 



n Hc 1-23 x 10- 25 / He ,4472 



™H 6 



6.16 x 10- 26 / Hfi 



= 0.15 



(4) 



This corresponds to a He mass fraction, Y=0.35 (assum- 
ing a metal fraction of 0.02). Such high He fractions 
are proposed in the second generation stars of globular 
clusters. However, we caution that the density and tem- 
perature of the gas that is assumed in this calculation 
are not very well constrained. While the result is rela- 
tively insensitive to higher temperatures at this density, 
at higher densities and temperatures the predicted line 
ratios are consistent with primordial abundances. 

3.2.3. The nature and future of Till 

The stellar populations in the two Till clusters are 
consistent with the multiple populations seen in the 
Galactic globular clusters: they have similar metallici- 
ties; an age difference of 70 Myr; and there is even the 
suggestion that the younger cluster is He enhanced. As 
demonstrated above, their proximity and similar veloc- 
ities are consistent with the clusters being bound. The 
coalescence of clusters like Tlll-1 and Tlll-2 may play 
an important role in the history of at least some star 
clusters. The resulting mixed population could produce 
a globular like cluster with multiple populations, similar 
to those observed in old globular clusters. Mergers of 
cluster systems hav e been propose d in environments like 
the Antennae (e.g. iKroupal 119981) . Additionally, merg- 
ers have been proposed as a possible explanation for ob- 
servations of some star clusters - such as the extended 
and massive Milky Way halo globular cluster NGC 2419 
(Brvins & Kroupa 2011), the 'faint fuzz y' star clusters 
observed in the halos of other galaxies (jBriins et al.ll2009l : 
lAssmann et al.ll2011[) and the obse rved mass segregation 
in some Galactic open clu sters (e.g. lMcMillan et alj2007t 
IMoeckel fc Bonnelll 120091 ). 

The coalescence of sub-clusters clusters may be impor- 
tant in the history of at least some clusters. However, 
the Till system is clearly different than the prototypi- 
cal multiple population cluster that is proposed by most 
current models. These models generally consider isolated 
clusters where the second populations form in the center 
of the clusters (though in reality the situation is likely 
more complex). Identification of such multiple popula- 
tions, within an individual cluster, may be possible in 
spectroscopic surveys of more clusters in the age range 
of lOs-lOOs of Myr. 

4. CONCLUSIONS 

We have considered the integrated emission from 
star clusters composed of single and multiple stel- 
lar populations. The stellar population models of 
iBruzual fc Cha riot (2003) identify two periods in a clus- 
ter's lifetime when the optical U-B and V-I colors of a 



multiple stellar population are significantly different from 
the colors of a SSP. However, we note that observational 
uncertainties, reddening and uncertainties in the metal- 
licity of a cluster make it difficult to confirm the presence 
of multiple populations in a cluster based on these two 
optical colors alone. Additionally, we note that there are 
significant variations in the colors predicted by different 
stellar population models at young ages. This further 
complicates the analysis of cluster colors. We note that 
better understanding of variations between models and 
the extension of this work to the ultraviolet and near- 
infrared may help to identify multiple population clusters 
based on broadband colors in the future. 

We also find that it is challenging to identify multiple 
stellar populations (with small age differences of order 
100 Myr) in the spectra of old clusters. This is because 
the spectral features of the different populations tend to 
be similar and their optical spectra are adequately rep- 
resented by SSP models. While future work may poten- 
tially be able to identify multiple populations based on 
inconsistent line ratios for a SSP, we note that such stud- 
ies will require well understood stellar population models 
and high signal to noise spectra. 

We have identified a strong and clear signal of multiple 
stellar populations in a young cluster's spectrum. This 
occurs when the youngest stellar population is <10 Myr 
old. The young massive stars present at these young ages 
will produce strong nebular emission. Such emission can 
not be produced by star clusters with ages <16 Myr. 
Therefore, a cluster featuring both strong absorption 
lines (that require an older stellar population) and nar- 
row nebular emission lines (that require a young popu- 
lation of stars) could only be produced by two distinct 
stellar populations. 

As a preliminary search, we have looked for such a clus- 
ter in published spectroscopy of 15 massive young star 
clusters in the Antennae. Seven of these clusters show 
emission lines and are consistent with young stellar pop- 
ulations and seven show absorption lines that are con- 
sistent with older ages. Interestingly, the other cluster 
shows both strong absorption and emission lines (sugges- 
tive of two populations with ages of 8 Myr and 80 Myr). 
We show that this emission comes from two neighboring 
sub-clusters that are consistent with being bound. The 
properties of these two clusters suggest that their coales- 
cence could produce a single multiple population cluster 
that is similar to that of the Galactic globular clusters. 

If all massive star clusters contain multiple stellar pop- 
ulations, as appears to be the case, then spectroscopy of 
more clusters in the 40-200 Myr range should identify a 
cluster with its second population younger than 10 Myr 
and hence showing these dual spectral features consistent 
with multiple stellar populations. Identification of these 
clusters is clearly extremely desirable as it can help to 
constrain the still poorly understood formationary the- 
ory of star clusters in the light of multiple stellar popu- 
lations. 
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